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摘 要 :了 ERAS 是 由 欧洲 中 期 天 气 预 报 中 心 (European Centre for Medium-Range Weather Forecasts, 
ECWMF) 推 出 的 新 一 代 再 分 析 产 品 , 可 为 地 面 观测 站 点 稀少 地 区 提供 新 的 降水 数据 来 源 。 基 于 
2008 一 2017 年 内 蒙古 45 个 地 面 站 点 逐日 降水 数据 ,采用 多 项 评价 指标 评估 了 ERA5 再 分 析 降 水 数 
据 在 多 时 空 尺 度 上 的 精度 ,并 结合 13 个 极端 降水 指数 和 综合 加 权 模 型 ,采用 Sen 斜 率 法 和 Mann- 
Kendall 趋势 检验 法 分 析 了 1981—2021 年 该 地 区 极端 降水 及 其 危险 性 时 空 变 化 规律 。 结 果 表 明 : 
(1) ERAS 能 够 较 好 地 再 现 降水 过 程 ,与 观测 数据 在 月 . 季 尺 度 上 呈 极 强 相 关 性 (相关 系数 CC> 
0.85), 在 日 尺度 上 表现 为 强 相 关 (CC=0.68), 在 空间 尺度 上 ,ERA5 数据 在 东部 站 点 探测 精度 优 于 西 
部 。(2) 除湿 日 降水 强度 (SDII) 、 强 降水 总 量 (R95pTOT) 以 及 持续 干燥 日 数 (CDD) 外 ,其 余 极 端 降 
水 指数 在 内 蒙古 均 呈 现下 降 趋 势 , 湿 日 总 降水 量 (PRCPTOT) 下 降 趋 势 最 快 ,为 -15.74 mm (10a)'。 
(3) 极端 降水 指数 在 空间 上 呈现 明显 区 域 分 异 规律 ,极端 降水 在 内 蒙古 西部 表现 为 强度 增加 ,在 中 
部 表现 为 频次 、 强 度 和 持续 时 间 均 减少 ,在 东部 则 表现 为 强度 增加 、 频 次 和 持续 时 间 减 少 。(4) 极端 
降水 危险 性 指数 在 鄂尔多斯 市 .呼伦贝尔 市 . 巴 放 淖尔 市 和 兴安 盟 等 人 口 相对 密集 、 经 济 发 展 较 快 


的 市 ( 盟 ) 存 在 高 值 中 心 , 呈 显著 上 升 趋势 ,应 给 


a D sA 
重点 关注 。 
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随 着 全 球 气候 变 暖 ,提高 降水 数据 对 极端 降水 
事件 的 监测 精度 ,是 提升 气 修 变化 研究 水 文 模拟 
预测 能 力 以 及 灾害 风险 管理 水 平 的 关键 环节 。 内 
蒙古 自治 区 (简称 内 莹 证) 地 处 干旱 和 半 干 旱季 风 
气候 区 ,地 幅 了 辽阔, 但 地 面 气象 站 点 分 布 稀少 且 不 
均匀 ,使 其 降水 时 空 特征 研究 受到 很 大 限制 。 
此 ,找到 适合 内 蒙古 的 高 分 辨 率 网 格 化 降水 资料 对 
本 区 域 开展 气候 变化 降水 特征 研究 及 灾害 风险 管 
理 具 有 重要 意义 。 近 些 年 来 ,许多 学 者 尝试 基于 地 
面 气象 站 点 和 遥感 卫星 降水 产品 等 数据 对 极端 降 
水 开展 了 一 系列 研究 。 其 中 ,气象 站 点 资料 最 为 准 
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确 , 但 观测 资料 时 间 长 度 有 限 , 且 无 法 反映 出 降水 
的 连续 空间 分 布 特征 ;还 感 卫星 降水 产品 包括 热 
带 降雨 测量 任务 TRMM'*" 全 球 降水 观测 GPM” 
气象 预测 中 心 校正 后 的 产品 CMORPH”、 基 于 人 工 
神经 网 络 估算 的 降水 产品 PERSIANN-CDR" 等 ,该 
类 数据 能 较 直 接 反 映 降 水 的 时 空 特征 ,可 解决 地 面 
站 点 数据 缺失 问题 ,但 卫星 反 演 得 到 的 降水 数据 原 
理 与 地 面 实测 不 同 ,需要 地 面 观测 数据 进行 验证 和 
改良 才能 更 好 使 用 ”; 再 分 析 资 料 , 如 ERA5S 降水 数 
据 等 ,综合 了 地 面 站 点 观测 .遥感 反 演 方法 的 优点 ， 
是 近 几 年 被 广泛 使 用 的 数据 类 型 。 


基金 项 目 : 国家 社会 科学 基金 重大 项 目 (18ZDA105 ) ;华东 师范 大 学 地 理 信息 科学 教育 部 重点 实验 室 主 任 基 金 (KLGIS2022C01) ; 国家 自 


然 科 学 基金 项 目 (42001222 ) ; 冰冻 圈 科 学 国家 重点 实验 室 开 放 基 金 项 目 (SKLCS 2020-02 ) 资 助 


作者 简介 : 牛 怡 莹 (2000-) , 女 ,硕士 研究 生 ,主要 从 事 多 灾害 综合 风险 评估 和 适应 研究 . E-mail: 1193418540@qq.com 
通讯 作者 : 王 军 (1975-) , 男 ,博士 ,教授 ,主要 从 事 城市 灾害 风险 与 公共 安全 管理 研究 . E-mail: jwang@geo.ecnu.edu.cn 


202309.00231v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


9 期 牛 怡 莹 等 : 内 蒙古 ERA5 再 分 析 降 水 数据 性 能 评 佑 与 极端 降水 时 空 特征 分 析 1419 


欧洲 第 五 代 ERA5 再 分 析 资 料 是 采用 最 先进 的 
全 球 资 料 同 化 系统 和 数据 库 , 对 各 种 来 源 的 观测 资 
料 ( 卫 星 地面 观 测 飞机 船舶 .无线 电 探 空 及 测 风 
气球 等 ) 进 行 质量 控制 和 同化 处 理 得 到 的 ,具有 时 
间 序 列 长 .分辩 率 高 等 优点 二 。 目 前 ,已 有 研究 对 
比 评 估 ERA5 再 分 析 降 水 数据 与 其 他 降水 数据 集 的 
性 能 :如 Chen 等 利用 4 种 再 分 析 降 水 产品 揭示 了 
东亚 地 区 暧 季 的 日 变 率 ,并 证 实 再 分 析 产 品 可 用 来 
描述 降水 的 日 循环 ; Gabriela 等 中 对 比 评 估 了 
ERAS ERA-Interim 以 及 CFSR 等 再 分 析 数 据 在 美 
南部 和 墨西哥 地 区 对 降水 指标 的 模拟 性 能 ,发现 
ERAS 明显 优 于 其 他 2 组 数据 ; 温 婷 婷 等 “对 比 评估 
T CRU, ERA5,CMFD 4 种 降水 数据 在 青藏 高 原 的 
适用 性 ,得 出 ERA5 数 据 与 观测 值 呈现 极 强 相关 性 ， 
年 总 降水 量 相关 系数 超过 0.9。 上 述 研 究 均 证 明了 
ERAS 再 分 析 降 水 数据 相 较 其 他 降水 数据 性 能 
优 ,能 够 准确 重 现 降 水 时 空 分 布 ,但 同时 也 发 现 
ERAS 存 在 高 估 降 水 强度 的 问题 55。 

AS SCX ERAS 再 分 析 降 水 数据 与 内 蒙古 地 面 站 
点 观测 降水 数据 进行 比较 ,从 多 时 空 尺度 评价 
ERA5 在 内 蒙古 的 适用 性 ;基于 此 ,参考 世界 气象 组 
织 推 荐 的 13 个 极端 降水 指数 ” ,探讨 研究 区 1981 一 
2021 年 极端 降水 的 时 空 变 化 规律 ,并 对 研究 区 极端 
降水 危险 性 时 空 变化 特征 进行 分 析 。 本 研究 可 为 
内 蒙古 等 地 面 监 测 站 点 稀少 地 区 气候 变化 降水 和 


N 图 例 

À @ 地 面 气象 站 点 
一 一 国界 线 
[| ] 省 界 
C] 盟 市 边界 


高 程 /m 


灾害 风险 评估 等 研究 提供 数据 来 源 ,同时 可 为 该 区 
域 制 定 气 候 变 化 适应 和 防汛 减灾 策略 提供 科学 
依据 。 


1 研究 区 概况 


Dy S2 Ti ( 37°24’ ~53°23'N , 97°12’ ~126°04' E) Ht 
处 中 国 北部 边疆 ,面积 118.3x10' km ,地 势 由 东北 
向 西南 延伸 , 横 跨 东北 华北、 西北 地 区 , 共 包 含 呼 
和 浩特 市 、 鄂 尔 多 斯 市 等 在 内 的 12 个 地 级 市 
( 盟 )。 气 候 以 温带 大 陆 性 气候 为 主 ,年 平均 降水 量 
50~450 mm ” , 且 年 内 分 布 不 均 , 由 西南 向 东北 逐渐 
增加 ,春季 多 大 风 天 气 ,夏季 短促 炎热 , 易 发 生 极端 
降水 ,秋冬 季 多 发 生 霜 冻 .寒潮 天 气 。 内 蒙古 地 势 
较 高 ,平均 海拔 1000 m 左 右 ( 图 1)。 


2 数据 与 方法 


1 数据 来 源 

地 面 台 站 观测 资料 为 国家 气象 信息 中 心 提 供 
的 中 国 地 面 气候 资料 日 值 降 水 数据 (http://data.cema. 
cn/) ,选择 2008 年 1 月 1 日 一 2017 年 12 月 31 日 作为 
评估 时 段 ,经 严格 的 质量 控制 , 共 筛 选 出 内 蒙古 45 
个 地 面 站 点 。 

ERAS 再 分 析 降 水 数据 源 于 欧洲 中 期 天 气 预报 
中 心 (ECMWF) 第 五 代 全 球 再 分 析 资 料 ERA5 数据 


注 :该 图 基于 国家 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 审 图 号 为 GS(2022)4312 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 
1 内 蒙古 气象 站 点 分 布 


Fig. 1 Distribution of meteorological stations in Inner Mongolia 


202309.00231v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1420 tid wA 46% 


集 (https://eds.climate.copernicus.eu/ ) 。 人 研究 选用 的 
数据 起 止 时 间 为 1979 年 1 月 1 日 一 2021 年 12 月 31 
日 ,时 间 分 辨 率 为 逐 小 时 ,空间 分 辩 率 为 0.23" ,日 、 
月 斥 度 数据 分 别 由 次 一 级 斥 度 累加 得 到 ,日 太 度 降 
水 数据 事先 已 进行 时 间 校 正 。 

2.2 研究 方法 


间 的 线性 相关 程度 ,最 优 值 为 1; 相 对 偶 差 (BIAS ) 用 
来 衡量 数据 的 系统 偏差 程度 ,最 优 值 为 0; 均 方 根 误 
差 (RMSE) 用 于 衡量 数据 之 间 的 离散 程度 ,最 优 值 
为 0。 分 类 指标 中 ,命中 率 (POD) 用 来 表现 待 评估 
数据 集 准 确 捕捉 实际 降水 事件 的 能 力 ,最 优 值 为 1; 
误 报 率 (FAR ) 用 来 反映 竺 评估 数据 集 的 误 报 情况 ， 


2.2.1 ERAS 再 分 析 降 水 数据 评估 方法 “降水 数据 
的 评估 方法 主要 有 2 种 :一 种 是 “ 格 点 ”对 “ 格 点 ”的 
评估 ,即将 网 格 数据 插值 到 地 面 站 点 上 ,与 地 面 站 
点 实测 值 进行 比较 ”*; 男 一 种 是 “ 格 网 ”对 “ 格 网 ” 
的 评 佑 ,往往 通过 特定 的 空间 化 方法 比如 克 里 金 插 
值 双 线性 插值 等 ”” ,将 站 点 数据 插值 得 到 面 状 数 
据 再 进行 比较 。 本 文采 用 双 线 性 插值 方法 ,考虑 站 
点 周围 的 4 个 网 格 作为 ERA5 在 该 站 点 上 的 数据 ,分 
别 在 年 . 季 月 .日 4 种 时 间 尺 度 上 ,对 ERA5S 与 地 面 
站 点 降水 数据 进行 “ 格 点 ”对 “ 格 点 ”的 评估 。 其 中 ， 
季 尺 度 的 评估 依据 我 国 气象 部 门 传统 的 划分 方法 ， 
将 4 季 划 分 为 春季 (3 一 5 月 ) .夏季 (6 一 8 A) 秋季 
(9 一 11 月 ) 冬季 (12 月 一 翌年 2 月 )。 

评估 指标 选取 方面 ,选用 了 误差 评价 指标 和 分 
类 统计 指标 共计 6 个 指标 ”( 表 1)。 在 评估 前 , 首 
先 对 ERA5 再 分 析 降 水 数据 进行 时 间 校 正 。 统 计 指 
标 中 ,相关 系数 (CC) 代 表 实 测 数据 与 待 评估 数据 之 


表 1 评价 指标 描述 


Tab. 1 Description of evaluation indicators 


评价 指标 计算 公式 最 优 什 
>C-OG -5 
相关 系数 (CC) = ~ 1 
[Sc x Is-5 
S,-G, 
相对 偏差 (BIAS)/% = x 100% 0 
$c 
均 方 根 误差 (RMSE )/mm 0 
命中 率 (POD ) H+M 1 
误 报 率 (FAR) HF 0 
关键 成 功 指数 (CSD) p 1 


注 : n 为 样本 数量 ; $ 为 ERA5 再 分 析 数 据 ; G, 为 地 面 站 点 观测 数 
据 ; S$、G 分 别 为 ERA5 和 地 面 站 点 观测 数据 的 均值 ; 万 为 ERA5 准 
确 捕捉 到 降水 事件 数 ;M 为 ERA5 漏 报 的 降水 事件 数 ;为 ERA5 误 
报 的 降水 事件 数 。 


最 优 值 为 0; 关 键 成 功 指 数 (CST) 综 合 POD 和 FAR 情 
况 , 反 映 待 评估 数据 集 监 测 实际 降水 事件 的 真实 能 
力 ,最 优 值 为 1。 
2.2.2 极端 降水 指数 计算 方法 ”本 研究 从 极端 降水 
的 频率 、 强 度 和 持续 时 间 3 个 维度 ,参考 世界 气象 组 
织 定 义 的 13 个 极端 降水 指数 来 表征 降水 的 极端 态 
过 程 ” ,用 于 分 析 内 蒙古 极端 降水 事件 的 时 空 变 化 
特征 ,所 选 指标 定义 见 表 2。 

另外 ,本 研究 综合 13 个 极端 降水 指数 ,构建 了 
极端 降水 危险 性 指数 (R) ,以 此 来 评价 内 蒙古 暴雨 
洪涝 灾害 危险 性 31。 


R= YX, (1) 
i=1 


式 中 : w, 为 第 i 个 极端 降水 指数 的 权重 ;X, 为 第 i 个 
极端 降水 指数 ;i 取 1,2,3,…, 13。 指 标 对 应 权重 采 
用 主 成 分 分 析 法 确定 ( 表 2)。 

2.2.3 ae EMER MSH Theil-Sen Median 方 
法 又 称 为 Sen 斜率 估计 ,是 一 种 稳健 的 非 参数 统计 
的 趋势 计算 方法 。 该 方法 对 于 测量 误差 和 离 
群 数 据 不 敏感 ,适用 于 长 时 间 序 列 数据 的 趋势 分 
析 , 斜率 为 正 表 示 呈 增加 趋势 ,斜率 为 负 则 表示 呈 
减少 趋势 。Mann-Kendall 趋势 检验 (简称 M-K 检 
验 ) 可 以 有 效 检 验 气象 水 文 要 素 的 长 时 间 序 列 的 趋 
势 变 化 和 显著 性 ,和 常用 于 降水 和 干旱 频次 等 趋势 的 
检测 。 本 文 取 显 著 性 水 平 为 0.05 , 若 检 验 量 Z> 
0, 则 表示 序列 旦 上升 趋 势 , 且 |2|>|Zu|= 1.96, 表 
示 序 列 呈 显著 上 升 趋势 ; 同 理 , 奉 检验 量 Z<0, 是 
|Z] >| Zo] = 1.96, 表示 序列 呈 显 著 下 降 趋 势 ,UF 和 
UB 统 计量 交点 表示 序列 的 时 间 突 变 点 。 


3 结果 与 分 析 


3.1 ERAS 降 水 数据 精度 评价 
3.1.1 不 同时 间 尺 度 精 度 评 估 

(1) 年 尺度 

基于 2008 一 2017 年 ERA5 降 水 数据 和 同期 内 蒙 
古 45 个 地 面 站 点 降水 数据 反 距离 权重 插值 结果 ,得 
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表 2 极端 降水 指数 定义 
Tab.2 Definition of extreme precipitation indices 
指数 指标 中 文 名 称 定义 单位 权重 
频率 指数 R10 中 雨天 数 一 年 内 日 降水 量 =>10 mm 的 天 数 d 0.0690 
R20 暴雨 天 数 一 年 内 日 降水 量 =20 mm 的 天 数 d 0.0845 
R30 特大 暴雨 日 数 一 年 内 日 降水 量 =30 mm 的 天 数 d 0.0868 
SDII 湿 日 降水 强度 年 总 降水 量 与 日 降水 量 =1 mm 天 数 的 比值 mm: d 0.0897 
强度 指数 R95p 强 降水 量 一 年 内 日 降水 量 > 第 95% 分 位 值 的 降水 量 之 和 mm 0.0909 
R99p 极 强 降水 量 一 年 内 日 降水 量 > 第 99% 分 位 值 的 降水 量 之 和 mm 0.0962 
PRCPTOT 湿 日 总 降水 量 一 年 内 日 降水 量 >1 mm 的 降水 量 总 和 mm 0.0700 
R95pTOT 强 降水 总 量 R95p 与 PRCPTOT 的 比值 % 0.0916 
RxI day 1 d 最 大 降水 量 一 年 内 最 大 日 降水 量 mm 0.0950 
RxSday 5 d 最 大 总 降水 量 一 年 内 连续 5 d 最 大 累积 降水 量 mm 0.0919 
Rx7day 7 d 最 大 总 降水 量 一 年 内 连续 7 d 最 大 累积 降水 量 mm 0.0904 
持续 时 间 CWD 持续 湿润 日 数 一 年 内 日 降水 量 =1 mm 的 最 长 持续 天 数 d 0.0340 
HA CDD 持续 干燥 日 数 一 年 内 日 降水 量 <1 mm 的 最 长 持续 天 数 d 0.0100 


到 内 蒙古 多 年 平均 年 降水 量 空间 分 布 图 (图 2):2 种 
数据 得 到 的 降水 量 等 级 分 界线 基本 重合 ,降水 量 空 
间 分 布 特征 一 致 , 均 反应 出 内 蒙古 降水 由 西南 向 东 
北 递增 的 特点 ,年 平均 降水 量 均 在 50~600 mm; 呼 伦 


N 
(a) 气象 站 点 人 


图 例 降水 量 /mm (| | 
— BRA C] ER CO 盟 市 边界 e690 500 400 300 200 100 20 


贝尔 市 降水 最 充沛 ,年 平均 降水 量 多 达 550 mm; 阿 
拉 善 盟 西 部 降水 最 少 , 年 平均 降水 量 不 足 100 mmo 


另外 ,内 蒙古 东部 赤峰 市 ERA5 的 降水 量 相 较 于 地 
面 站 点 降水 量 存在 低 佑 。 


2 气象 站 点 与 ERA5 年 平均 降水 量 空间 分 布 


Fig.2 Spatial distributions of annual average precipitation of meteorological stations and ERAS 


(2) ERE 

由 ERA5 降 水 量 与 地 面 站 点 降水 量 的 季 尺 度 散 
点 密度 图 (图 3a~d) 可 知 :ERA5 数 据 在 夏季 表现 最 
好 ,CC 达 0.89,BIAS 为 0.02,RMSE 为 59.36 mm; i 
在 冬季 表现 最 差 ,CC 为 0.82,BIAS 为 0.72,RMSE 为 
8.59 mm, RMSE 的 差异 与 夏季 降水 量 多 而 冬季 降水 
少 有 关 。 另 外 从 BIAS 指标 可 以 发 现 ,ERA5S 降水 数 
据 在 内 蒙古 存在 高 估 ,这 可 能 与 ERA5 数 据 本 身 算 
法 误差 有 关 。 


(3) 月 日 尺度 

ERAS 与 地 面 站 点 逐 月 降水 数据 变化 趋势 一 
致 ,ERA5 降水 数据 在 夏季 存在 一 定 误 差 , 对 最 大 月 
总 降水 量 存在 一 定 程度 高 估 ( 图 4) ,这 与 月 尺度 散 
点 密度 图 (图 3e) 结 果 一 致 ,但 在 2013 年 和 2014 年 8 
月 ERA5 低 估 了 月 降水 量 。ERA5 降水 数据 与 站 点 
数据 在 日 尺度 呈 强 相关 性 (CC=0.68) , 低 于 月 尺度 
(CC=0.88), 且 在 2 种 时 间 尺 度 上 ERA5S 均 存在 高 估 
(图 3f)。 
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图 3 不 同时 间 尺 度 的 ERA5 
Fig. 3 Density-colored scatterplots of ERAS and precip 
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图 4 ERA5 与 气象 站 点 降水 数据 的 月 均 降 水 量 对 比 


Fig.4 Comparison of monthly average precipitation between ERAS and precipitation data of meteorological stations 


31.2 基于 站 点 尺度 的 精度 评估 基于 站 点 尺度 对 
ERAS 降水 数据 进行 精度 评价 ,以 0.1 mm: d“ BEAK VE 
为 是 否 发 生 降水 的 冰 值 交 ,结果 如 图 5 所 示 。 从 误 
差 评 价 指标 看 ,CC 在 东部 地 区 站 点 表现 优 于 西部 ， 
在 日 尺度 上 48.9% 的 站 点 CC 为 0.70~0.80; 线 性 相关 


性 大 的 站 点 往往 也 表现 出 BIAS 高 的 特点 , 除 5 个 地 
面 站 点 表现 为 较 小 程度 的 低估 以 外 ,其 余 站 点 BIAS 
均 表 现 出 大 于 10% 的 高 估 ; 降 水 量 大 的 地 区 RMSE 
更 大 ,RMSE 呈现 出 由 西南 向 东北 递增 的 趋势 ,在 数 
值 上 由 1~2 mm #4 £ 3~5 mm。 从 分 类 统计 指标 看 ， 
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(a) CC À (b) POD À 


CC 
@ <0.50 @ 0.60~0.70 
@ 0.50~0.60 ©® 0.70~0.80 


(c) BIAS A 


BIAS 
@ <0.00 


@ 0.10~0.20 


(e) RMSE 


RMSE/mm 
@ 4.00~5.00 @ 2.00~3.00 


POD 
@ 0.60~0.70 @ 0.80~0.90 
@ 0.70~0.80 @>0.90 


(d) FAR A 


FAR 
@ >0.50 @ 0.40~0.45 
@ 0.45~0.50 @ <0.40 


(Ð CSI 


图 例 ”一 一 国界 线 Lo) 省 界 LCL) 盟 市 边界 


注 :CC 为 相关 系数 ;POD 为 命中 率 ;BIAS 为 相对 偏差 ;FAR 为 误 报 率 ;RMSE 为 均 方 根 误差 ;CSI 为 关键 成 功 指数 。 下 同 。 
图 5 ERA5 在 不 同 误差 评价 指标 与 分 类 统计 指标 下 探测 能 力 表 现 


Fig. 5 Detection capability performance of ERAS under different error evaluation indicators and classification statistics indicators 


POD 在 内 蒙古 东部 表现 优 于 西部 ,西部 站 点 能 捕捉 
到 90% 以 上 的 降水 事件 ,东部 站 点 POD 为 0.80~ 
0.90 ,整个 内 蒙古 仅 有 3 个 站 点 的 POD 在 0.70 之 下 ; 
75.6% 的 站 点 FAR 为 0.40~0.50, 且 西部 站 点 大 于 东 
部 站 点 ,这 与 东西 部 降水 量 分 布 不 均 有 很 大 关系 ; 
CSI 反 应 出 ERA5 数据 集 监测 实际 降水 事件 的 真实 
能 力 , 与 POD 指 标 表 现 一 致 。 

3.1.3 不 同 强度 降水 的 探测 精度 评估 为 了 评估 
ERAS 降水 数据 对 于 不 同 强 度 降水 探测 精度 的 差异 
程度 ,绘制 了 POD、FAR 和 CSI3 个 分 类 统计 指标 在 
不 同 国 值 下 的 折线 图 (网 6)。POD AI CSI MA BE 
的 增 大 而 减 小 , 即 ERA5 降水 数据 的 探测 性 能 随 降 
水 阔 值 的 增加 而 下 降 ,FAR 整 体 呈 小 幅 上 升 的 态 
势 。0.1 mm-d :党 被 作为 判断 是 否 发 生 降 水 的 阔 


值 ,此 时 POD 可 达到 0.9, 即 ERA5 能 捕捉 到 90% 的 


降水 事件 ;而 在 阅 值 达到 15 mm dE, ERAS 精度 
明显 下 降 ,FAR 会 高 于 POD。 


0.8 — POD — FAR -« CSI 


01 1 3 5 8 10 12 15 18 20 22 
降水 阔 值 /mm*d-1 


6 WEKE EE F ERAS 降水 数据 探测 能 力 表 现 
Fig. 6 Detection capability performance of ERAS 


precipitation data under different precipitation thresholds 
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3.2 极端 降水 指数 时 空 演变 特征 

3.2.1 极端 降水 指数 的 室 间 变化 13 个 极端 降水 指 
数 在 内 蒙古 均 呈 现 出 自 西向 东 的 空间 分 布 特征 (图 
7) ,其 中 频率 指数 R10、R20 和 R30 自 西 向 东 依 次 递 
减 , 呼 伦 贝 尔 市 降水 量 东 西 差 异 最 大 , 横 跨 5 个 等 
级 ;从 降水 强度 来 看 ,呼和浩特 市 .通辽 市 .呼伦贝尔 
市 等 SDI 可 达到 7~8 mm*d-', R95pTOT 可 达 24%~ 
28%; 持 续 时 间 指 数 CDD 在 25~200 d 范 围 内 ,东西 跨 
度 大 , 除 阿 拉 善 盟 外 ,其 余地 区 CDD 均 少 于 100 d, 
CWD 空间 分 布 特征 与 CDD 指数 相反 , 且 在 数值 上 


人 


(a) R10 
R10/d 国 8~12 


m0~4 E 12~16 
m48 EE 10~20 


(d) SDN 人 (e) R95p 


SDII/mm-d"! m 5~6 
mm 2~4 Ml 6~7 
mi 4~5 mi 7~8 


(g) PRCPTOT A (h) R95pTOT 


PRCPTOT/mm i 300~400 
ma 0~150 m 400~500 


mm 150~300 mm 500~680 E 16~20 


N 
(j) RxSday A (k) Rx7day 


Rx5day/mm m 40~60 
EE 10~20 mi60~80 


m 20~40 E 80~100 30~350 


(m) CDD 
CDD/d m80~100 


R95p/mm m 70~100 
m4~30 m100~130 
30~70 E 130~180 


R95pTOT/% = 20~22 
ml0~16 =m 22~24 
E 24~28 


Rx7day/mm = 50~70 
ml10~30 im 70~90 
EE 90~110 
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wi 4 46% 


远 小 于 CDD。 

从 极端 降水 空间 变化 趋势 来 看 (图 8) ,频率 指 
数 R10 在 内 蒙古 东部 呼伦贝尔 市 .赤峰 市 .通辽 市 、 
锡林郭勒 盟 以 及 呼和浩特 市 等 市 ( 盟 ) 呈 明显 下 降 
趋势 ,最 大 Sen 斜率 为 -1.667 mm (10aj ,表现 出 在 
人 口 密集 区 域 明显 降低 的 特征 ,R20 仅 在 呼和浩特 
市 有 明显 下 降 趋势 ,而 R30 在 整个 内 蒙古 变化 趋势 
不 显著 。 强 度 指 数 R95p、Rx5day、Rx7day 在 内 蒙古 
西部 阿拉 善 盟 、 鄂 尔 多 斯 市 以 及 东部 呼伦贝尔 市 局 
部 地 区 旦 增加 趋势 ,在 中 部 呼和浩特 市 .锡林郭勒 


人 


(c) R30 

R30/d 

m0- 2-3 
ml~2 3-4 


A (f) R99p A 


R99p/mm m 25~35 
mel~10 mm 35~45 
mm10~25 mm 45~65 


A (i) Rxlday 


Rxlday/mm m 30~40 
7~20 40~50 
20~30 m50~60 


N 
人 (1) CWD 


CWD/d m 5~7 
mm3 7-9 
m35 E 9~12 


mE 25~60 m 100~150 图 例 

60~80 国 150~200 — 国界 线 
[| 省 界 
C 最 市 边界 


注 : 图 中 各 指标 缩写 词 含 义 见 表 2。 下 同 。 
图 7 1981 一 2021 年 13 个 极端 降水 指数 的 空间 分 布 
Fig.7 Spatial distributions of 13 extreme precipitation indices from 1981 to 2021 
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图 8 1981 一 2021 年 13 个 极端 降水 指数 的 空间 变化 趋势 
Fig. 8 Spatial variation trends of 13 extreme precipitation indices from 1981 to 2021 


盟 .赤峰 市 以 及 通辽 市 等 部 分 地 区 呈 减 少 趋势 ,其 
HH R95p HY Sen 斜率 最 大 ,为 19.47 mm-(10a)'; SDI, 
R95pTOT 以 及 Rxlday 在 内 蒙古 西部 阿拉 善 盟 、 鄂 尔 
多 斯 市 .中 部 锡林郭勒 盟 以 及 东部 赤峰 市 .通辽 市 、 
兴安 盟 和 呼伦贝尔 市 等 局 部 地 区 呈 增 加 的 趋势 ,而 
在 呼和浩特 市 .赤峰 市 等 局 部 地 区 呈 减 少 趋 势 ; 
PRCPTOT 在 内 蒙古 大 部 分 区 域 均 呈现 显著 减 小 趋 
势 ;R99p 在 内 蒙古 变化 趋势 不 显著 。 持 续 时 间 指 数 
CDD 在 内 蒙古 西部 阿拉 善 盟 ,中 部 巴 彦 淖尔 市 ` 包 
头 市 . 乌 兰 察 布 市 ,以 及 东部 呼伦贝尔 市 下 降 趋势 
显著 ;CWD 在 锡林郭勒 盟 东部 和 呼伦贝尔 市 北部 地 


区 下 降 趋势 显著 。 

3.2.2 极端 降水 指数 的 年 际 变化 ”1981 一 2021 年 极 
端 降水 的 年 际 变化 的 频率 .强度 及 持续 时 间 变 化 趋 
势 如 图 9 所 示 , 除 CDD 外 ,其 余 指数 在 内 蒙古 均 表 
现 为 下 降 趋势 。 频 率 指数 R10、R20、R30 在 近 40 a 
的 Sen 斜率 分 别 为 -0.439 d- (10a)! ,-0.125 d:(10a)” 
和 -0.036 d.(10aj,R30 变 化 趋势 小 的 原因 为 内 蒙古 
日 降水 量 难以 达到 30 mm。 强度 指 数 ROSpTOT 呈 
增加 趋势 ,Sen 斜率 为 0.32%.(10aj 5;SDI 呈 微弱 增 
加 趋势 ,Sen 斜率 为 0.001 mm -+ d-' + (10a)'; R95p、 
R99p、PRCPTOT、Rxlday、Rx5day 以 及 Rx7day J) Œ 
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9 1981 一 2021 年 13 个 极端 降水 指数 的 年 际 变化 趋势 和 M-K 检 验 结 


Fig.9 Inter-annual variation trends and M-K test results of 13 extreme precipitation indices from 1981 to 2021 


减少 趋势 ,但 下 降 趋 势 不 显 著 , 其 中 PRCPTOT 下 降 
速率 最 快 ,为 -17.36 mm.(10aj-:。 频 率 和 强度 指数 
均 在 1999 年 存在 时 间 突 变 点 ,1999 年 前 后 趋势 变化 
程度 或 方向 发 生 明 显 变 化 ,可 见 1999 年 是 内 蒙古 气 
候 变 化 的 关键 时 间 点 。 持 续 时 间 指 数 CWD 和 CDD 
空间 变化 相反 ,其 Sen 斜率 分 别 为 -0.150 d- (10a) 和 
4.924 d.(10aj, 且 2 个 指数 在 2010 年 均 有 时 间 突 变 
点 ,在 2010 年 前 后 时 间 序 列 呈 不 同 程度 的 增加 ( 减 
小 ) 趋 势 。 
3.3 极端 降水 危险 性 时 空 演变 特征 

由 1981 一 2021 年 内 蒙古 极端 降水 危险 性 指数 的 
时 空 变化 特征 (图 10) 可 以 发 现 , 内 蒙古 极端 降水 危险 


H REKEDT . 乌 海 市 .阿拉 善 盟 等 地 极端 降水 危 
分 性 指数 呈 上 升 趋势 , 且 在 阿拉 善 盟 东部 存在 高 值 
中 心 ;而 在 呼和浩特 市 .马兰 察 布 市 .包头 市 . 锡 林 
郭 勒 盟 、 赤 峰 市 .通辽 市 等 地 呈 下 降 趋势 , 并 在 马兰 
察 布 市 存在 低 值 中 心 。 


4 讨论 


可 靠 的 降水 估算 不 仅 是 气候 变化 影响 趋势 和 
波动 等 研究 的 重要 内 容 , 也 是 开展 干旱 、 洪 涝 等 灾 
害 监 测 与 应 对 的 重要 基础 性 工作 。 从 论文 采取 的 6 
个 统计 指标 来 看 ,ERA5 再 分 析 降 水 数据 在 内 蒙古 


性 指数 整体 呈 下 降 趋 势 ,Sen 斜率 为 -0.041.(10aj 1， 
极端 降水 危险 性 指数 在 1999 年 发 生 显著 突变 ， 
1981—1999 年 和 2000 一 2021 年 危险 性 呈 上 升 趋 
势 。 内 蒙古 极端 降水 危险 性 指数 在 空间 上 存在 显 
著 差 异 : 有 56% 区 域 极端 降水 危险 性 指数 呈 上 升 趋 
势 , 其 中 1.93% 的 区 域 呈 显著 上 升 , 极 端 降水 危险 性 


具有 和 较 高 的 可 靠 性 ,可 作为 内 蒙古 降水 研究 的 一 种 
新 的 数据 源 。 该 数据 在 日 尺度 上 的 相关 性 最 低 , 这 
可 能 是 因为 年 . 季 和 月 降水 量 是 由 日 降水 量 合成 ， 
存在 误差 抵消 效应 。 这 种 效应 与 东亚 地 区 ”青藏 
高 原 地 区 “以 及 长 江 三 角 洲 地 区 ”等 相关 人 研究 具有 
一 致 性 。 研 究 也 发 现 EFRAS5 存 在 高 估 问 题 ,推测 可 能 
是 因为 部 分 降水 在 被 站 点 检测 到 之 前 被 蒸发 ,也 


指数 的 Sen 斜 率 为 -0.157~0.177:(10a)'!; 在 内 蒙古 东 
北部 的 呼伦贝尔 市 和 兴安 盟 , 西 南部 的 巴 彦 淖尔 


可 能 与 再 分 析 技 术 中 格 网 对 流 参数 化 方案 的 缺陷 
和 对 固体 降水 的 不 准确 估计 有 关 ,但 总 体 上 误差 
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Fig. 10 Spatiotemporal variation characteristics of extreme precipitation danger index in Inner Mongolia from 1981 to 2021 


2 

极端 降水 的 变化 会 诱发 城市 时 小 灾害 ,威胁 城 
市 安全 ,影响 社会 经 济 发 展 。 本 研究 在 前 人 研究 基 
础 上 引入 了 极端 降水 危险 性 指数 ,将 极端 降水 与 其 
危险 性 建立 联系 ,揭示 了 内 蒙古 极端 降水 危险 性 的 
时 空 变化 特征 ,研究 结果 与 白 美 兰 ”、 尤 莉 等 “和 刘 
泓 志 等 ”的 研究 结果 相 一 致 。 尽 管内 蒙古 极端 降水 
的 频次 整体 呈 减 少 趋势 ,但 在 局 部 地 区 极端 降水 的 
强度 却 出 现 明显 增加 趋势 。 内 蒙古 极端 降水 危险 


5 结论 


二 


(1) ERAS 降水 数据 能 较 好 再 现 内 蒙古 降水 过 
程 , 但 相 较 于 站 点 观测 数据 ,整体 存在 高 估 。 月 尺 
度 (CC=0.88) 精 度 表 现 优 于 日 尺度 (CC=0.68) ,这 是 
由 于 月 尺度 削减 了 日 尺度 上 的 降水 差异 ; 季 尺 度 
E , ERAS 与 站 点 观测 降水 数据 在 夏季 相关 性 最 高 ， 
而 在 冬季 最 小 ,这 可 能 与 夏季 降水 量 大 .冬季 降水 


性 在 人 口 相 对 密集 、 经 济 较 发 达 的 市 ( 盟 ) ,如 哪 尔 
多 斯 市 .呼伦贝尔 市 . 巴 彦 浑 尔 市 、 兴 安 盟 表现 为 增 
加 趋势 ,应 关注 极端 降水 对 这 些 城市 带 来 的 洪 游 灾 
害 影响 ,提升 极端 降水 预报 预警 和 洪涝 灾害 应 急 响 
应 能 力 。 另 外 ,部 分 区 域 极端 降水 危险 性 指数 显著 
下 降 , 如 锡林郭勒 盟 草 原 等 ,这 一 结果 与 金 令 等 ”、 
马 梓 策 ”等 研究 得 出 的 锡林郭勒 盟 草 原 东部 属于 
旱灾 高 危险 性 的 认识 相 一 致 。 内 蒙古 草原 地 区 降 
水 趋 于 减少 可 能 会 对 草原 生物 多 样 性 和 生态 屏障 
功能 产生 影响 ,进而 影响 到 区 域 生态 安全 , 故 应 关 
注 这 些 区 域 气候 变化 趋势 预 估 ,并 加 强 降水 监测 与 
生态 风险 预警 。 


量 小 有 关 。 从 空间 上 看 ,ERA5 数据 在 东部 站 点 的 
表现 明显 优 于 西部 , 且 ERA5 的 探测 性 能 随 着 降水 
强度 的 增加 而 下 降 。 

(2) 从 极端 降水 指数 的 年 际 变化 来 看 , 近 40 a 
内 蒙古 年 湿 日 总 降水 量 (PRCPTOT) 减 少 ,Sen 斜率 
为 17.36 mm (10a)', KR oh BE FE BL SDI, ROSpTOT 以 
及 持续 时 间 指 数 CDD 呈 上 升 趋势 ,其 余 极 端 降水 指 
数 R10, R20, R30, R95p, R99p, Rxlday, Rx5day, 
Rx7day、CWD 均 呈现 下 降 趋势 ; 除 CDD 外 ,其 余 12 
个 极端 降水 指数 均 在 1999 年 发 生 时 间 突 变 ,1999 年 
为 内 蒙古 气候 状态 发 生变 化 的 关键 点 。 

(3) 从 空间 上 看 ,内 蒙古 极端 降水 呈现 明显 的 
“西南 一 东北 " 带 状 分 布 ,表现 为 由 西南 向 东北 递 
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增 , 空 间 变 化 趋势 表现 出 明显 的 区 域 分 异 特征 。 内 
蒙古 西部 地 区 除 部 分 强度 指数 (SDIT、PRCPTOT、 
Rx5day、Rx7day) 呈 现 增加 外 ,其 余 指 数 均 无 明显 变 
化 趋势 ;中 部 地 区 多 数 极端 降水 指数 (R10、R20、 
SDII, R95p, PRCPTOT,, R9SPTOT Rxlday、 Rx5day、 
Rx7day 以 及 CDD) 呈 显著 减少 趋势 ;东部 地 区 则 在 
多 数 强 度 指数 (SDIT、R95p、R99p、R95PTOT、Rxlday、 
Rx5day 以 及 Rx7day) 上 表现 为 增加 ,而 在 频率 指数 
(R10) 和 持续 时 间 指 数 (CWD CDD ) 上 表现 为 减少 。 

(4) 近 40 a 内 蒙古 极端 降水 的 频次 整体 呈现 减 
少 趋势 ,但 在 局 部 地 区 极端 降水 的 强度 却 呈 现 显著 
增加 趋势 。 综 合 13 个 极端 降水 指数 发 现 ,极端 降水 
危险 性 指数 在 驾 尔 多 斯 市 .呼伦贝尔 市 ` 巴 彦 济 尔 
市 、 兴 安 盟 等 人 口 相对 密集 、 经 济 发 展 较 快 的 市 
( 盟 ) 存 在 高 值 中 心 ,表现 为 显著 增加 趋势 ,Sen 斜率 
H 0.177- (10a)! ,未 来 应 关注 这 些 区 域 的 极端 降水 变 
化 特征 ,加强 监测 .预警 与 应 急 响应 。 
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Abstract: ERAS is a new generation reanalysis product launched by the European Center for Medium- Range 
Weather Forecasts that can provide a new source of precipitation data for areas with few ground observation sta- 
tions. Based on the daily precipitation data of 45 ground stations in Inner Mongolia, China from 2008 to 2017, we 
evaluated the accuracy of ERAS reanalysis precipitation data on multiple temporal and spatial scales using multi- 
ple evaluation indicators and applied a comprehensive weighting model to construct an extreme precipitation dan- 
ger index by integrating 13 extreme precipitation indices. We then used principal component analysis, Sen’ s 
slope method, and a Mann-Kendall trend test to analyze the temporal and spatial changes of extreme precipitation 
and extreme precipitation danger in the region from 1981 to 2021. The results show that: (1) ERAS can reproduce 
the precipitation process better; the ERAS precipitation reanalysis dataset performs better at a monthly time scale 
than at a daily time scale, and its accuracy is highest in summer and lowest in winter. ERAS precipitation data are 
highly correlated with observed data at monthly and seasonal scales (correlation coefficient>0.85) and are strong- 
ly correlated with observed data at a daily scale (correlation coefficient=0.68). The detection accuracy of ERAS 
data is better for the eastern stations than for the western stations. (2) Extreme precipitation indices showed a 
downward trend except for daily wet precipitation intensity, total heavy precipitation, and continuous dry days. 
The annual total wet day precipitation declined fastest, and the Sen’ s slope value was — 15.74 mm (10a) '. (3) 
The extreme precipitation indices show clear regional spatial differentiation. Extreme precipitation shows an in- 
crease in intensity in western Inner Mongolia, a decrease in frequency, intensity, and duration in mid-Inner Mon- 
golia, and an increase in intensity, frequency, and duration in eastern Inner Mongolia. (4) The extreme precipita- 
tion danger index has high central values and shows significant upward trends in Ordos City, Hulun Buir City, 
Bayannur City, Hinggan League, and other cities/leagues with relatively dense populations and rapid economic 
development and should therefore receive special attention. The results of this study can facilitate the discovery 
of better datasets for the analysis of climate factors in Inner Mongolia, and the study provides a theoretical basis 
for the formulation of climate change adaptation measures and future climate prediction. 


Key words: ERAS; precipitation; accuracy evaluation; extreme precipitation index; Inner Mongolia 


